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ABSTRACT 

Red-giant stars are proving to be an incredible source of information for testing models of stellar 
evolution, as asteroseismology has opened up a window into their interiors. Such insights are a 
direct result of the unprecedented data from space missions CoRoT and Kepler as well as recent 
theoretical advances. Eclipsing binaries are also fundamental astrophysical objects, and when coupled 
with asteroseismology, binaries provide two independent methods to obtain masses and radii and 
exciting opportunities to develop highly constrained stellar models. The possibility of discovering 
pulsating red giants in eclipsing binary systems is therefore an important goal that could potentially 
offer very robust characterization of these systems. Until recently, only one case has been discovered 
with Kepler. We cross-correlate the detected red-giant and eclipsing-binary catalogs from Kepler 
data to find possible candidate systems. Light-curve modeling and mean properties measured from 
asteroseismology are combined to yield specific measurements of periods, masses, radii, temperatures, 
eclipse timing variations, core rotation rates, and red-giant evolutionary state. After using three 
different techniques to eliminate false positives, out of the 70 systems common to the red-giant and 
eclipsing-binary catalogs we find 13 strong candidates (12 previously unknown) to be eclipsing binaries, 
one to be an non-eclipsing binary with tidally induced oscillations, and 10 more to be hierarchical triple 
systems, all of which include a pulsating red giant. The systems span a range of orbital eccentricities, 
periods, and spectral types F, G, K, and M for the companion of the red giant. One case even suggests 
an eclipsing binary composed of two red-giant stars and another of a red giant with a J-Scuti star. 
The discovery of multiple pulsating red giants in eclipsing binaries provides an exciting test bed for 
precise astrophysical modeling, and follow-up spectroscopic observations of many of the candidate 
systems are encouraged. The resulting highly constrained stellar parameters will allow, for example, 
the exploration of how binary tidal interactions affect pulsations when compared to the single-star 
case. 

Keywords: stars: AGB and post-AGB — (stars:) binaries: eclipsing — stars: oscillations (including 
pulsations) — stars: interiors — techniques: photometric — methods: data analysis 



1. INTRODUCTION 

Binary stars are prime targets to study stellar evolu- 
tion: for some spectral types more than 50% of stars are 
estim ated to belong to binary or multiple-sta r systems 



(e.g., Duquennoy & Mayor 1991 Lada 2006[), and Ke- 



pler's thTroTTiiwall^^ physical 
parameters, which is of crucial importance to constrain 
stellar evolution models. It is possible to determine the 
masses of each component (Mi, M 2 ) for visual or eclips- 
ing binaries, provided that spectral lines are detected for 
each star to track the Doppler shifts along their orbits, 
and the product (Mi + M2)sinz for spectroscopic-only 
binaries, where i is the orbital plane inclination. More- 
over, close-in binary systems host unique and poorly- 
understood physical processes such as mass exchange be- 
tween the two stars. 



The Kepler satellite ( |Borucki et al.|2010[ ) detected 2616 
eclipsing binaries (hereafter EBs) by the end of Jan- 



( here aft' 



uary 2013, which represente d about 1.4 % of th e targets. 
These ta rgets were listed by Prsa e t al. (20 11 ), then up - 



dated by Slawson et al. (2011 ) and Matijevic et al. (2012 ) 



Within this list, binary systems were classified into lour 
types, (i) Detached (D), where each component's radius 
is smaller than its Roche lobe so stars are spherical. The 
stars have no major effect on each other, and essentially 
evolve separately. Most binaries belong to this class, (ii) 
Semi-detached (SD), where the biggest star fills its Roche 
lobe leading to mass exchange. The mass transfer dom- 
inates the evolution of the system. In many cases, the 
inflowing gas forms an accretion disc around the accre- 
tor. (hi) Over-contact (OC), where both stars fill their 
Roche lobes and are in contact. The uppermost part of 
the stellar atmospheres forms a common envelope that 
surrounds both stars, (iv) Ellipsoidal variation (ELV), 
where no eclipse is observed but the system is detected 
in Kepler's light curves through the ellipsoidal shape of 
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the st ars. In addition to the Kepler list, Coug hlm et al.| 
(2011) proposed a list of low-mass (< 1 Mq) EBs, of 



which man y in fact belong to the lists of |Slawson et al. 
fl2011|) and|Matijevic et al.| ([2012]) . 



Red giants are evolved stars that have depleted the 
hydrogen in their cores and are no longer able to gen- 
erate energy from core hydrogen burning. The physical 
processes taking place in their interiors are fairly poorly 
understood. However, the study of the global pulsations 
of red giants (hereafter RGs) with asteroseismology is 
capable of unambiguously determining bulk properties 
such as mass, radius, temperature, metallicity, and also 
the evolutionary state of RGs. Indeed, t he new era of 
space-based missions such as CoRoT (Bagl in et al.| 2QQ9) 
and Kepler has dramatically increased the amount and 
quality of the available asteroseismic data. In particu- 
lar, global oscillations of several thousands of solar-like 
stars and RGs have been detected in Kepler data, and 
analysis of the oscillation eigenmodes now allows robust 
seismic inf erences to be drawn about their internal struc- 
ture (e.g.JChaplin et all|2011| |Hekker et al,||2QQ9| |Bed-| 
ding et al.l|2(JlU| IMosser et al.||20T0| |Huber et aTp OlO ; 
Bedd ing et al.||201ip . 

Applying modern asteroseismology to EB systems for 
which photometric and radial velocity data exist leads 
to the best possible physical characterizations: this is 
because masses and radii may be measured in two inde- 
pendent ways. Such systems are cornerstones for testing 
stellar evolution models. Until recently, only one case of 
an oscillating RG belonging t o an EB syst em had been 
reported (KIC 8410637, pekker et al.|2010[ ). As only one 
eclipse of that system was observed, no estimate of or- 
bital period or eccentricity could be obtained, but the 
global oscilla tions of the RG sta r were clearly detected. 
In addition, |Derekas et al.| ( |2Q11| ) report the detection of 
a triple system containing a KG (HD 181068), and they 
explicitly mention that solar-like oscillations are not vis- 
ible, even though most stars with similar parameters in 
the Kepler database do clearly show such oscillations. 
In general, the Kepler mission has succeeded in making 
breakthroughs in both the fields of binary stars and as- 
teroseismology. For example, Kepler helped reveal the 
presence of tidally-induced pulsations in the binary sys- 
tem KOI 54, which are the result of resonances between 
the dynamic tides at periastron a nd the free oscill ation 
modes of one or both of the stars (W elsh et al.|201l] ). 

Since most stars are observed at a 29-min cadence with 
Kepler, global modes of main-sequence solar-like stars 
are not accessible; however, global modes of RG stars 



larger than 3.5 Rq are accessible (e.g., Table 3 of Mosser 



et al. 2012c). Fortunately, an RG catalog from the Kepler 



learn has been compiled and made public for the scien- 
tific community. 1 In this paper, we establish a list of 
RG candidates that likely belong to EBs or multiple-star 
systems, which we obtain from the EB and RG public 
catalogs. We test whether these candidates are part of 
EB or multiple systems and characterize their main phys- 
ical properties. Note we do not work with radial velocity 
measurements, for which data acquisition is in progress. 
We first present the results of a cross-correlation of the 
EB and RG catalogs (Section [2|, and the subsequent de- 
tailed analysis of their light curves to determine eclipse 



and asteroseismic properties (Section [3|. We identify 70 
RGs possibly belonging to EB systems, of which 47 show 
clear global oscillation modes. Mean properties of the 
global modes are used to infer RG masses and radii. We 
study several ways to determine whether the oscillating 
RGs actually do belong to the EBs with which they are 
associated, and then describe details of several impor- 
tant cases in Section |4j In Section [5] we conclude by 
defining the observations that are needed to fully char- 
acterize this set of stars and discuss general implications 
for oscillating stars in binary systems. 

2. DATA 
2.1. Kepler data 

All data used in this paper are photometric measure- 
ments obtained by the NASA satellite Kepler, launched 
in M arch 2009 to search f or exoplanets in the habitable 
zone ( Borucki et al.||2010 ). Since its launch, Kepler has 
been monitoring about 156 000 stars in a 105 deg 2 field 
of view in between the Cygnus and Lyra constellations. 
Data are subdivided in quarters Q, i.e., three- month runs 
at the end of which the satellite is rotated by 90 deg to 
maintain the Sun's position on its solar arrays and to 
keep the radiator pointed to deep space. The commis- 
sioning quarter Q0 and quarter Ql lasted respectively 
10 and 35 days, respectively. Here, we only utilize public 
data that are available through Q13 as of January 2013. 

Light curve data released for the public on the MAST 
database can be obtained in raw or corrected form. Cor- 
rected fluxes have been processed by the Presearch Data 
Conditioning (PDC) pipeline that removes signatures in 
the light curves that are correlated with systematic error 
sources from the telescope and spacecraft, such as point- 
ing drift, focus changes, and thermal transients. The 
PDC attempts to correct for such errors while preserv- 
ing planet transits and other astrophysically interesting 
signals. Further details of the pipeline are desc ri bed by 



Kinemuchi et al.| ([2012]), |Stumpe et al.| ([2012]), |Smith 
et al.| ( |2012[ ), and the Kepler handbook. 2 We use cor- 
rected PDC fluxes in this work. 

In addition, for each object there are target pixel files, 
consisting of the flux for all the pixels contained within a 
predefined mask which are used to create the data found 
in the photometric light curve files. Each target pixel 
file contains these pixels as a time series of images in a 
bin ary FITS table. We use target pixel files in Section 
|3.1| to determine whether the eclipsing signal we observe 
is from a nearby contaminating object. 

The Kepler observations are sampled at either long 
cadence (29.4244 minutes) or short cadence (58.89 sec- 
onds, for only 512 targets). Most targets in this study 
were observed at long cadence; for the few objects that 
were also observed at short cadence it was never for more 
than 30 days. Therefore, we work primarily with long 
cadence data, and only consider the short-cadence data 
to search for high-frequency modes of the main-sequence 
or subgiant companion star belonging to the considered 
EB system. For long-cadence data, the Nyquist cut-off 
frequency is ^ nyq = 283.2 /iHz, which limits the possibil- 
ity of performing asteroseismology on any object whose 



frequency at maximum amplitude v n 



By as- 



1 http://archive.stsci.edu/kepler/red_giant_release.html 



2 http:/ /archive. stsci.edu/kepler/manuals/archive_manual.pdf 
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suming an RG effective te mperature of 4800 K and l M fl , 
asteroseismic scaling laws (Kjeldsen fc Bedding 1995 Hu- 
ber et al. 2011 Belkacem et al. 2011) predict that the 
lower limiting radius is approximately 3.5 Rq. 

2.2. Red-giant and binary-system catalogs 

The Kepler team released a list of 13 698 RGs on 
September 27, 2011, selected from the target list us- 
ing color-magnitude estimates and considering effective 
temperatures T e ^ = 4800 ± 300 K and surface gravities 
logg = 2.9 ± 0.5 (cgs). This latter criterion cuts off RGs 
with predicted oscillation peak frequency larger than 
about 320 /iHz, slightly above the long-cadence Nyquist 
frequency. We note that not all of the RGs were contin- 
uously observed, and they are typically the first targets 
to be dropped from the survey when pixel resources be- 
come scarce (see the Kepler red giant database). The 
RG Kepler magnitudes range from 7.9 to 14.0 which 
makes them sufficiently bright for ground-based spec- 
trometry to determine atmospheric parameters (as the 
SDSS III APOGEE/ Kepler experiment APOKASC at 
A pache Point Obse rvatory is currently realizing) . 



Prsa et al. ( 20111) released a catalog of EBs identified 
in the Kepler field from the first data releases (Q0-Q1). 
The catalog was motivated by the exquisite quality of 
Kepler data which has led to the discovery of hundreds 
of new systems, revolutionized accuracy in modeling EB 
systems, and provided an estimate of the frequency of oc- 
currence of EBs. Their method uses the Transit Planet 
Search (TPS) algorithm, first developed to detect ex - 



2010a|bD , 



oplanets in the Kepler data (Jenkins et al 
by adapting it to search for eclipse durations consistent 
with EBs. The targets that present a positive detec- 
tion are then down-selected to exclude objects already 
identified as exoplanets, variable or spotted stars that 
mimic eclipse shapes, blends from background stars, or 
pointing jitter artifacts. The eclipses are subsequently 
modeled with EBAI (Eclip sing Binaries via Artificial In- 
telligence; |Pri^r^r][2008| , which relies on trained neu- 
ral networks using synthetic eclipse profiles to extract 
the physical parameters. For contact systems, orbits are 
assumed to be circular and the main physical parame- 
ters fit are the temperature ratio T2/X1, the mass ratio 
q — M2/M1, the fillout factor F which is a function of 
the potential gravity, and the inclination sinz. For semi- 
detached and detached binaries, the mass ratio cannot be 
estimated, and F is replaced with the sum of fractional 
radii p\ + p2 = (Ri + R2)/a 1 where a is the semi-major 
axis of the binary orbit. The eccentricity e together with 
the argument of the periastron uj are further introduced 

through e cos uo and e sin uj. 

A second catalog from |Slawson et al.| ( |2011[ ) updates 
the results by including data from Q2, adding some Ke- 
pler Objects of Interest (KOIs) flagged as possible exo- 
planets and later determined to be EBs, adding objects 
with period longer than 44 days (i.e., Q0+Q1), rejecting 
variables stars initially considered to be EBs, and remov- 
ing EBs that were initially blends situated at the edge 
of the photometric aperture later re-observed with a re- 
centered aperture. This updated catalog presents 2165 
EBs, composed of 58% detached (D), 7% semi-detached 
(SD), 22% overcontact (OC), 6% ellipsoidal (ELV), and 
7 % undetermined systems. These studies conclude that 
eclipsing binaries represent 1.4% of the Kepler target 




Figure 1. From left to right and top to bottom, histograms of Ke- 
pler magnitude, contamination factor, effective temperature, and 
surface gravity for the entire EB catalog (dashed line, left y-axis) 
and for the sample flagged as RG and EB (solid line, right y-axis). 
Data are from the KIC. 

list. A third update has been published that e mploys an 
automatic m orphological classification scheme (Matijevic 



et al^2012 ), and additional EB candidates have been re 
leased on the Kepler MAST EB database, leading to a 
total of 2616 EBs. 

2.3. Cross- correlation of both lists 

The Kepler database contains 13 698 stars identified 
as RGs, while the eclipsing binary catalog contains 2616 
targets, as of January 2013. The cross-correlation of both 
catalogs reveals that 70 stars are flagged as both EB and 
RG (hereafter RG/EB): 56 systems (66%) are classified 
as D, seven (10 %) as SD, 11 (16 %) as OC, and two (4 %) 
as ELV; two (3 %) are unclassified. Such proportions are 
close to those of the whole sample, which is initially sur- 
prising since we do not expect RGs to belong to close-in 
systems, in particular SD and OC, where orbital periods 
range from 0.2 to 6 days (see TablepT). For our purpose of 
considering stars flagged as EB ana RG, it is prudent t o 
keep in mind the recommendations of |Prsa et al. (2011): 



1. "The high star density leads to a non- negligible 
likelihood of associating an EB event with the 
wrong star;" 

2. "The EB interpretation should be taken with ex- 
treme caution for stars with high fraction of flux 
contamination;" 

3. "Stars with shallow eclipse events should also be re- 
garded with caution even if the flux contamination 
is modest." 

Hereafter, when we mention RG/EBs, we imply candi- 
date RG/EB systems. To understand the specificities of 
RG/EBs with respect to the whole EB catalog and how 
likely they may correspond with misidentifications, we 
compare histograms of observing conditions and stellar 
atmospheric parameters (estimated when systems were 
suppo sedly single stars) f rom the Kepler Input Catalog 
(KIC," 
lished 



Brown et al. 201 lL and orbital parameters pub- 
by |Slawson et al. "j[2011 ) in Figures Jl] and |5J 



It appears that RG/EBs are slightly brighter on aver- 
age than the full EB catalog, but this is not significant 
since RGs were selected with a magnitude limit of 14. 
For stellar crowding, we only show contamination factors 
calculated on the apertures corresponding to the second 
of the four positions of the satellite (i.e., during Ql, Q5, 
and Q9). We also tested the three other satellite orien- 
tations without noting any relevant difference. We see 
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Figure 2. From left to right and top to bottom, histograms of or- 
bital period, sum of relative radii, stellar temperature ratio, orbital 
plane inclination sini, and eccentricity for the whole EB catalog 
(dashed line, left y-axis) and for the sample flagged as RG and EB 
(solid line, right y-axis). The x-axis is truncated at an eccentricity 
of 0.6 because so few sys tems are above this value. Data are from 
the |Slawson et al.] ( |2011| ) EB catalog. 

that the contamination factor of RG/EBs is lower than 
the mean contamination of all EBs, which shows that our 
candidates are not overly predisposed to the risk of tar- 
get misidentification or blending. Histograms of RG/EB 
effective temperature and surface gravity are not repre- 
sentative of the whole EB catalog, but are instead typical 
of RGs. This means that even if a RG does not belong 
to an associated EB, there is still one RG per photo- 
metric aperture, whose flux overwhelms that from any 
secondary star. We note two exceptions at temperatures 
around 6500 K (see Table [T]) that are classified as OC 
and are likely misidentified as RGs. 

The orbital period histograms of the RG/EBs and the 
total EB sample in Figure [2] are similar, and shows that 
most RG/EB systems have orbital periods < 10 days. 
The sum of relative radii are also similar in both sam- 
ples with the exception of the lack of values lower than 
0.1 in RG/EBs. The relative deficit of stars with a tem- 
perature ratio equal to one among RG/EBs can be re- 
lated to the minor proportion of contact systems in the 
RG/EB sample. We note that the inclination and or- 
bital eccentricity distributions seem to be locked in the 
ranges s'mi = [0.99,1] and e = [0.1,0.15], which suggests 
that t he estimate of these parameters by Slaw son et al.| 
(2011) could be biased. Measuring orbital eccentricity 
from photometric data alone, i.e., without radial veloci- 
ties, is known to be inaccurate. 

To get an initial rough estimate whether a RG is com- 
patible with an EB, we determine the amplitude of the 
deepest eclipse measured in the RG/EB light curve. We 
assume the system is observed from the orbital plane, 
and that the RG is larger than its companion. We there- 
fore adopt the convention that the deepest eclipse is the 
secondary eclipse (i.e., companion star going behind the 
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Figure 3. Expected primary eclipse depth (given by the contours, 
in %) by assuming a companion star of temperature T2 and radius 
R2 eclipsing a RG of temperature 4800 K and radius 10 Rq. 

RG) if the companion's temperature is higher than the 
RG's, and that the deepest eclipse is the primary eclipse 
otherwise (i.e., when the RG's temperature is higher than 
that of the companion star). Primary and secondary 
eclipses will have equal depths when the two stars' effec- 
tive temperatures are equal. A raw proxy of relative pho- 
tometric dimming during primary and secondary eclipses 
can be obtained by neglecting stellar limb darkening and 
comparing simple luminosities: 



SI 



^2 1 2 



primary 



H 2 1 2 



1 J secondary n l 1 1 "•" n 2 1 2 



(i) 



where the subscripts 1 and 2 stand for the RG and the 
companion star, respectively. Note the numerator of the 
second equation is R%Tf since since the light dimming is 
due to the star of radius R2 which hides a surface ttR^ of 
the star of temperature Ti . In Figure [3] we present the ex- 
pected eclipse depth assuming a typical RG with radius 
10 Rq and effective temperature 4800 K. We conclude 
that the range of photometric dimmings ~ 0.2 — 20%, 
which we measure in all but nine of the RG/EBs, is 
compatible with a stellar pair that includes a RG and 
a main-sequence star (see Table [T]). We note that pri- 
mary eclipse depths will be lower than the theoretical 
proxy in Equation ([!]) when eclipses are grazing due to 
limb darkening. For the nine systems with very shallow 
eclipses, their depths range from 0.02 % to 0.14 %. These 
cases could either correspond to RGs with a background 
EB, EB systems with grazing eclipses, or EB systems 
where the companion star's size is a few percent of the 
RG's size, as could be possible for brown dwarfs or giant 
planets. 

From the RG and EB catalogs, as well as from the asso- 
ciated KIC parameters log g and T e ff , it is assumed that 
each candidate system contains one EB and one RG in 
the photometric aperture. However, suspicion of blend- 
ing arises from quick calculations regarding the orbital 
period distribution. Let us consider a simple configura- 
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Figure 4. Orbital velocity contours (km s — 1 ) of a RG in a bi- 
nary system as a function of the orbital separation and orbital pe- 
riod. A circular orbit is assumed. We further assume (Ri,M±) = 
(1(XR©,1.4M©) for the RG and (R 2 ,M 2 ) = (LR©,1M©) for the 
companion s tar. Red crosses repr esent the values of the RG/EB 
sample from Slawson et al. (2 011| ). The thick black line indicates 
the critical velocity tor a phase-locked RG. 

tion where a RG with radius and mass (10i?©, 1.4M ) 
is in a phase-locked binary system with a solar-like com- 
panion (li?©, IMq) in a circular orbit. (Such parameters 
for the RG are close to the median values of R G/EBs ob- 
tained with asteroseismology in Section [3~3j ) In Figure 
|4j we show the expected orbital velocities that the RG 
would have in that configuration, by varying both orbital 
period P and semi-major axis a: 



^orb,l 



2ira M 2 
~P~ Mi + M 2 " 



(2) 



By using the values for the sum of fractional radii and 
or bital periods est i mated for 46 out of the 70 systems 



by Slawson et al. (2011), we see that 14 RGs (mostly 
of the (J(J class) would have an orbital velocity between 
500 and 2600 kms" 1 . Even if our assumptions about 
the companion star's mass and radius are incorrect, this 
indicates that some systems are not physically possible 
as they would require an orbital velocity greater than 
the critical velocity that would begin to tear apart the 
RG. We consider a velocity to be "critical" when the 
rotation velocity at the equator equals the escape ve- 
locity. In the case of a phase-locked system, the rota- 
tion period is simply the orbital period. Any additional 
proper stellar rotation in the direction of the orbital mo- 
tion (e.g., Earth) would lower the critical velocity, while 
a stellar rotation in the opposite direction of the orbital 
motion (e.g., Venus, which is unlikely) would enhance 
the critical velocity. Based on this criterion, only 25 out 
of the 46 systems, for which we have orbital parame- 
ters, have a period compatible with the critical velocity 
threshold. This is the first of several clues that points to 
the fact that many of the candidates are not true RG/EB 
systems. 

In the following sections, we use several analysis tech- 
niques to better characterize the systems and more ac- 
curately determine how many of our candidate systems 
are bona-fide RG/EBs. 



3. LIGHT CURVE ANALYSIS FOR ECLIPSES AND 
ASTEROSEISMOLOGY 

3.1. Search for contamination from surrounding stars 

The cross correlation between the RG and EB 
databases resulted in 70 identifications. We consider four 
possible scenarios: 

1. The RG is actually one of the EB stars. 

2. The RG is aligned within the pixel field of view with 
the EB and is a part of a multiple system (gravita- 
tionally speaking), of which two close-in stars mu- 
tually eclipse, with the RG out of the EB's orbital 
plane. Eclipse Timing Variations (ETVs) may be 
observed in this case (see Section 4.1). 



3. The RG is aligned within the pixel field of view 
with the EB, but is not gravitationally bound. 

4. The RG and the EB fall on different pixels in the 
aperture and are not gravitationally bound. 

The last case can be verified using the target pixel files 
associated with each star by computing a map of the 
relative intensity variation dl/I to check whether the 
depth of the eclipse is correlated to the peak intensity 
source on the detector. 

Every quarter when the spacecraft rotates, a star falls 
onto a different set of pixels, and the pixel mask (aper- 
ture) that defines the optimal output light curve often 
changes. In the case of a long-period EB, the eclipses 
only occur in a few of the quarters. We detrended each 
pixel light curve in a quarter by a low-order polynomial 
and then modulated ( "folded" ) it by the period estimated 
from the global light curve. The relative intensity drop 
dl/I was then calculated and an image of this eclipse 
depth was compared to an average intensity image for 
the quarter. 

In 40 of the 47 cases, we find that the peak in intensity 
corresponds to the peak in eclipse depth, indicating sce- 
narios 1 or 3 above cannot be ruled out (and for the nine 
brightest stars, saturation in the aperture that buries any 
eclipse signature on individual pixels does not allow this 
method to be used). We find 7 clear cases where the pixel 
with the deepest eclipse occurs away from the pixel with 
the peak mean intensity (and typically outside of the Ke- 
pler aperture as well). An example is shown in Fig. [5] for 
KIC 4576968, which is classified in the binary catalog as 
an OC system. The overall peak in intensity of the mean 
photometry is well inside the defined aperture, yet the 
maximal eclipse depth is about two pixels outside of this 
region. The eclipse depth considering only pixels within 
the aperture is about 0.06%, while outlying pixels reach 
about 2%. Studying the individual folded light curves 
for each pixel confirms that the EB and the RG are on 
different pixels, and likely not part of the same system. 

All of these 7 cases occur in short, approximately less 
than 1 day period EBs, where it is indeed physically im- 
possible for the RG to be one of the eclipsing stars and 
still remain intact. Five such systems (KIC 2711123, 
4576968, 5652071, 7879404, 11968514) show behavior 
where we are confident that the RG is a background con- 
taminant and not actually a part of the EB, as scenario 
4 describes. Two other cases (KIC 7031714, 7955301) 
show ETVs and are consistent with scenario 2, but are 
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Figure 5. Identification of RG star KIC 4576968 and contam- 
inating EB from the Kepler Target Pixel Files. The top panels 
show (a) the map of the mean intensity, (b) the map of eclipse 
depth in each pixel, and (c) a segment of the light curve obtained 
by summing over all pixels (including outside the aperture). White 
corresponds to highest flux on a linear scale, and the red outlines 
the aperture used to compute the Kepler light curves. The ma- 
trix in panel (d) shows the period-folded light curve in each target 
pixel normalized to the median value. The gray boxes denote the 
aperture. The strongest eclipse signature is outside of the aper- 
ture, while the highest flux pixel is within the aperture. All data 
are from Q7. Note that the dispersion of the folded light curves is 
higher in the three pixels where the intensity is maximum because 
of the presence of RG oscillations. 

likely to be false positives. Indeed, if these objects were 
gravitationally bound and separated by one or more four- 
arcsecond pixels, the physical separation at the distance 
of a magnitude-9 giant star would be on the order of 
several hundred to thousands of AU, which implies an 
orbital period of several thousands of years. Hence these 
cases could correspond to a RG with a background triple 
system. 

In addition, for each star, we compute the power spec- 
trum of each pixel. It appears that the global oscillations 
always originate from the pixels where intensity is maxi- 
mum, which means that the brightest star systematically 
corresponds to a red giant, for the systems where oscil- 
lations are detected. We even distinguish by eye the RG 
oscillations on the map of folded light curves across the 
photometric aperture through a point dispersion larger 
than the average noise (see Figure |5j panel (d)). 

3.2. Cleaning the time series 

Two distinct ways of processing the light curves are 
required for modeling eclipses or analyzing seismic prop- 
erties. Eclipse modeling consists of extracting the stel- 
lar parameters from fitting the shape of the light curve 
during eclipses after removing instrumental biases and, 
ideally, other sources of stellar variability (e.g., activ- 
ity, pulsations). On the contrary, studying global mode 
properties requires that the time series from which we 
calculate power spectra is free from periodic signals such 
as instrumental systematics and stellar eclipses. In other 
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Figure 6. Cleaning the light curve of short-period binary KIC 
3532985. The top panel shows a segment of the original light curve 
(red) and the corrected light curve (black) after eclipse removal, 
which is used for seismology analysis. The middle panel shows the 
original light curve folded on the period (black) and its rebinned 
version (red). The bottom panel is the result after a smoothed 
light curve without eclipses is subtracted, and is used for light 
curve modeling. The dashed vertical lines in the top and bottom 
panel demarcate successive orbital periods. 

words: global modes are noise for eclipse modeling and 
vice versa. 

Nevertheless, there are initial reduction procedures 
common to both cases. For example, it is absolutely 
necessary to remove long time drifts and discontinuities, 
which are mostly of instrumental origin, particularly af- 
ter Kepler rotates between consecutive quarters. Note 
that we do not use the portion of each light curve between 
eclipses for our light curve models due to difficulties in 
disentangling instrumental flux variations, stellar activ- 
ity (spots, granulation), reflection, Doppler beaming, and 
ellipsoidal variations. We assume that the median fluxes 
of each quarter's data are equal, and normalize them to 
1 to work with relative fluxes. The entire light curve is 
then concatenated, and all outliers (selected as being out 
of the light curve mean dispersion by an amount evalu- 
ated individually for each target) are eliminated. The 
next step depends on the specifics of each light curve. 
Five classes of datasets and procedures are described: 

• Short orbital period EBs (e.g., Figure (6|. When 
the time series is long enough to contain more than 
about 30 orbits, we proceed by assuming that the 
signal fluctuations, be they stellar or instrumen- 
tal in origin, may be averaged out by folding and 
rebinning the light curve. We proceed by subtract- 
ing the mean folded light curve from the original. 
This results in a time series with no eclipse (see 
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Figure 7. Cleaning the light curve of long-period binary system 
KIC 8430105. The top panel is the same as Figure [6] and shows 
a segment of the original light curve (red) and the corrected light 
curve (black) after eclipse removal, which is used for seismology 
analysis. The middle panel shows a zoomed view of a single eclipse 
from the original light curve, and the polynomial fit used to re- 
move the eclipse signal for asteroseismology (red). The bottom 
panel is the result after a smoothed light curve without eclipses 
is subtracted, and is used for light curve modeling. The dashed 
vertical lines in the top and bottom panels demarcate the region 
highlighted in the middle panel. 

the black curve in Figure [6j top panel) which re- 
tains the stellar variability that we use to search 
for global modes in the Fourier domain. Next, 
we smooth this time series with a moving aver- 
age whose width is set on the characteristic time 
scale of the photometric variations to be cancelled. 
The smoothed time series is then subtracted from 
the original time series to get a light curve with 
a flat level between eclipses. In principle, this is 
the best method because it is simple and makes 
no assumption about the origin of any photometric 
fluctuations. In practice, we limited the use of this 
method to systems with orbital period shorter than 
20 days. 

• EBs with stellar variations on the order of the or- 
bital period, which is much longer than the eclipse 
duration (e.g., Figure |7|). In this case, a stellar 
signal of amplitude similar to or higher than the 
eclipse's depth cannot be cancelled out by fold- 
ing the data. We proceed by identifying the cen- 
ters of the primary and secondary eclipses, bridg- 
ing them with a second-order polynomial, and then 
subtracting the smoothed light curve from the orig- 
inal to yield a flattened light curve with equally 
deep eclipses. We use a second-order polynomial 
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Figure 8. Cleaning a light curve with variable eclipse depth and 
timing (KIC 7955301). The top panel shows the entire original 
light curve (red) and the corrected light curve (black) after eclipse 
removal, which is used for seismology analysis. The middle panel 
shows a zoomed view of one of the eclipses (black) and an example 
of the function used to fit each eclipse (red). The bottom panel 
shows the result after the eclipses are subtracted from the original 
light curve, and upon which asteroseismology is performed. No 
light curve modeling is attempted in these cases. 

to "fill" each eclipse by simply fitting the short re- 
gions of the light curve on either side of the eclipse 
of duration equal to the eclipses. We find a second 
order polynomial is sufficient account for the local 
behavior of the light curve. Next, an asteroseismic 
analysis is carried out on the light curve that has 
the eclipses filled with polynomials to avoid peri- 
odic gaps. This technique cannot be applied if the 
eclipse duration is too large compared to the or- 
bital period because we would cancel the signal by 
filling the eclipses. 

• Long orbital period EBs. When the eclipse du- 
ration is less than 5% of the duration of the en- 
tire light curve, we postulate that simply removing 
the data where eclipses occur does not significantly 
change the duty cycle for the asteroseismic analy- 
sis. For eclipse modeling in such systems, we flatten 
the light curve after "filling" eclipses with second- 
order polynomials, as in the previous case. 

• EBs with variable eclipse depth and timing (e.g., 
Figure [8]). In several cases, we encounter light 
curves with eclipses whose depth and timing change 
with time, either of astrophysical or instrumen- 
tal origin (e.g., varying flux contamination at each 
field rotation). We do not attempt to model such 
eclipses in this work because it often implies ac- 
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Figure 9. Cleaning an over-contact (OC) binary light curve with 
variable amplitude (KIC 7879404). The top panel shows the en- 
tire original light curve (black) and the smoothed version (red). 
Since the smoothing function is a weighted moving average, the 
smoothed light curve equals the signal at each interruption, giving 
an impression of spikes. The middle panel shows both the original 
light curve after subtraction by the smoothed version (black) and 
the local minima and maxima per half-period (red). The bottom 
panel is the result after the original light curve is divided by the 
smoothed maxima-to-minima distance. A zoomed view is is shown 
for time between 710 and 711 days. This final light curve is used 
for asteroseismology, and we attempt no light curve modeling. 

counting for a third body when there is an as- 
trophysical cause, or it requires careful modeling 
of systematic errors when the cause is instrumen- 
tal. For asteroseismology, however, we fit individ- 
ual eclipses with simple functions and then subtract 
them from the time series. The functions we use 
to fit eclipses are Gaussian when eclipses are graz- 
ing (no flat bo ttom, or "v-shaped"), or the |Mandel| 
& Agol ( 2002 ) exoplanetary fitting function other- 



• Contact EBs with variable amplitude (e.g., Fig- 
ure [9|. As indicated in the previous cases, varying 
flux contamination can modify the apparent ampli- 
tude of eclipses, and the technique of fitting eclipses 
with Gaussian or exoplanetary functions is not pos- 
sible for contact (OC) systems. Therefore, after 
subtracting a smooth version of the light curve, 
we measure positions and amplitudes of all local 
maxima and minima per half-orbit and smooth 
them. Then we normalize the time series with the 
smoothed amplitude of the light curve and sub- 
sequently measure asteroseismic parameters. No 
light-curve modeling is attempted for these sys- 
tems. 



3.3. Asteroseismic analysis 
3.3.1. Detection and properties of global pulsation modes 

We search for global oscillation modes in light curves 
whe re e clipses have been handled as described in Sec- 
tion 3.2 In this paper, we only consider clear oscillatory 
excess power, and do not consider low signal-to-noise ra- 
tio excess power that in principle could b e detected with 
filtered autocorr elation techniques (e.g., Mosser & Ap- 



pourchaux|[2009 ) . Indeed, since the light curves studied 
here could nave remnants of eclipse features, the presence 
of harmonics of the orbital period would strongly perturb 
the computation of the autocorrelation function, leading 
to false mode detections. 

Power density spectra (PDS) are computed with the 
discrete Fourier transform with no oversampling up to 
the Nyquist frequency (283 /iHz). We detect global 
modes in 47 of the 70 RG/EB candidates. The break- 
down by binary classification is 37 D, two SD, five OC, 
and two ELV. All power spectra showing solar-like oscil- 
lations are presented in Figure [To] 

One of the asteroseismic aims of this paper is the mea- 
surement of RG global mode parameters, such as their 
frequency at maximum amplitude ^ max and mean large 
frequency separation Av. The pair (^ max ? Ai/) allows 
us to estimate the RG's mass and radius through well- 
known asteroseism ic relations (Kjeldsen & Bedding 1995 
Huber et al.pOlO ). 



The frequency at maximum amplitude ^ max is obtained 
by fitting the PDS with a Gaussian function for the 
mode envelope and a sum of semi-Lorentzian functions to 
model the stellar activity. This standard technique was 
first suggested by 
widely used for so 



Harvey (1985) for the Sun, and is now 
ar-like stars (e.g., Chaplin et al.|2011 ). 



The PDS is divided by the semi-Lorentzian background 
to produce a whitened spectrum, and autocorrelation is 
then used to extract the mean large separation Aza Fig- 
ure [IT] shows Av versus ^ max for the RG/EBs compared 
witnTne empirical relationship that was established on 
thousands of red giants, subgiants, and main sequence 



stars with CoRoT and Kepler data (|Hekker et al.[ 2009 
Stell o et alT||2009l IMosser et al |[20iu| |Huber eTaflpHT T 
Mos ser et al.||2012a| ). Only the three stars with v n 
close to the Nyquist frequency are significantly different 
from expectations, because the determination of ^ max is 
biased by the PDS truncation. 

To compute RG masses and radii, estimates of their ef- 
fective temperatures T e ^ are required. As we co nclu ded 
from analyzing the T e ^ distribution in Section |2.3[ we 
can safely assume that temperatures from the Kepler 
database correspond the red giants for the RG/EB sam- 
ple. Using the asteroseismic parameters and these tem- 
per at ures^_aJong_w^^ asymptotic scaling factors 
for RGs (|Mosser et al.||2012b[), Figure 12 shows the esti- 



mated masses and radii lor the 47 pulsating RGs, which 
range from 0.8 — 2.8M and 3.7 — 33. Oi?©, respectively. 
The new scaling changes the values of mass and radius 
for each star by up to a few percent. These results are 
presented in Table [2] 

3.3.2. Mode identification and mixed modes 

Peaks in light curve power spectra correspond to 
unique stellar oscillations, and identifying these oscilla- 
tions yields enormous insight into interior stellar proper- 
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Figure 10. Power density spectra (PDS) of the 47 light curves in which global, solar-like oscillations are detected, sorted by increasing 
!/ m ax from top to bottom and left to right. Some spectra are truncated at low frequency because of the filtering applied to remove the 
signature of an eclipse. Note the different frequency x-axis scales in each column. Each system is labeled by its Kepler ID number. 
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Figure 11. Large frequency separation (Au) versus frequency at 
maximum amplitude (i/ max ) for the 47 pulsating RG/EB candi- 
dates derived from asteroseismology. The red line shows the em- 
pirical relationship e stablished from a large sample of RGs from 
|Mosser eta l. (2012a). The dashed vertical line denotes the Nyquist 
frequency tor the observation cadence, above which three RGs have 
their z/ max . 



^0 



I 



°20 



10 



+ ^ + 

++ /. + + 

++ 
+++ 




1.5 2 

Mass (M Q ) 

Figure 12. Histograms of masses and radii of the 47 pulsating 
RGs in RG/EB candidate systems. These have been derived astero- 
seismically a nd include new asymptotic scaling corrections ( |Mosser| 
|et al.||2012b| >. The scatter plot shows the correlation between the 
two observables. 
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Figure 13. Echelle diagrams of the 47 pulsating RG stars whose KIC names are labeled. The x-axis is the frequency modulo the large 
frequency spacing (i.e., from to Ai/), and the y-axis is frequency. The three noticeable dark mode ridges (for most stars) correspond to 
angular degree £ = 1, £ = 2, and £ = 0, from left to right, respectively indicated in dashed red, green and blue lines. Each case has been 
shifted so that the radial modes are roughly vertically aligned for illustration purposes. Each system is identified by its Kepler ID number. 



ties. For example, modes that are a mixture of acoustic 
and gravity modes provide diagnostics for probing the 
stellar core. Such mixed modes may be used to deter- 
mine if a star is on the red-giant branch (RGB) where 
it is still burning H in a shell surrounding the core, or if 
is part of the red clump (RC) after possibly experienc- 



ing a He flash and is now fusing He in its c ore (Bedding 
et al.|201Q||2Qll||Mosser et al.|2 Qll||2Q12c[). In addition, 
we may dist mguish the second 



may distinguish tne secondary red clump (RC2) that 
consists of stars more massive than 1.8 M©. These stars 
have started to burn He in a non-degenerate core, which 
disting uishes them fr om low-mass stars in the main red 
clump (Girardi 1999). These mixed modes typically ap- 
pear most clearly in the i — 1 modes as a "forest" of 
peaks. 

Figure [13] shows the echelle diagram of all 47 pulsating 
RGs in candidate RG/EBs. We detect the typical pair of 
ridges of solar-like oscillations (£ = 0, 2 and £ = 1 spaced 
by half the large separation) for 14 of the 47 oscillating 



RGs, while clear signatures of mixed £ — 1 modes are evi- 
dent in 23 stars. The remaining targets present confusing 
oscillatory spectra. These few RGs that show suppressed 
mixed modes certainly be long to the RG group identi- 



fied by Mosser et al. ( 2Q12a[ ). Such low- amplitude mixed 
modes nave been observed in a group of stars starting the 
ascent of the RGB, at all evolutionary stages. We note 
that 11 of the 14 stars with suppressed mixed modes have 
the longest orbital periods of all the RG/EBs. 

Careful analysis of the dipole mixed-mode signatures 
allows us to classify red giants as either RGB or RC 
stars. This was first co mputed by measuring the bumped 
mixed -mode spacings ( [Bedding et al.|2010 ; Mosser et al 
|2011|). However, using the asymptotic development of 



mixed modes provides the exact period spacin g AHi, 
which precis ely characterizes the radiative core (|Mosser| 



et al. 2012c). It also allows access to the core rotation 
rate ol these giants, and provides an estimate of P rot 



(Beck et al. 2012 Mosser et al. 2012b). The values of 
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these mixed- mode parameters for 12 of the 23 cases are 
given in Table [2] Interestingly, the cores are rotating 
with periods from 30 to about 385 days, seemingly un- 
corrected with their masses and increasing in most cases 
with radius. Also provided in Table [2] are the classifica- 
tions of 27 RGs. We note that when neither the gravity- 
mode period spacing nor the bumped spacing are avail- 
able, the value of the large separation Av can be used 
to classify RGs in certain cases: a RG with Av > 9 /iHz 
is on the RGB, but when 3.7 < Av < 5.5/iHz, the prob- 
ability of h aving a clump star is higher than 90% (see 
Figure 3 of |Mosser et al.|2012c[ ) 

3.4. Light curve modeling 
3.4.1. Modeling the light curves 

To independently compare physical parameters of a 
subset of the binary systems, we model their light curves 
using the Eclipsing Light Curve (ELC) code ( |Orosz| 
fc Hausch ildt 20 0£| an d/or the JKTEBOP code (e^7 
iSouthworth et al.|20U9 ). ELC in particular uses a genetic 
algorithm and Monte Carlo Markov Chain optimizers to 
simultaneously solve for a suite of stellar parameters. It 
is well-suited to this analysis because additional data can 
be added as it becomes available, such as radial velocities 
from spectra. Using only the Kepler light curves, we use 
ELC and JKTEBOP to solve for relative fractional radii 
Ri/a and R2/a, temperature ratio T 2 /Ti, orbital inclina- 
tion i, eccentricity e and e cos uj when applicable (where uo 
is the longitude of periastron), the Kepler contamination 
factor, and the stellar limb darkening parameters for the 
quadratic limb darkening law. The orbital period was de- 
termined previously through a simple iterative technique 
and held fixed during this analysis. In the limiting case 
where the two stars are sufficiently separated as to be 
spherical in shape, E LC has a " fast a nalytic" mode that 
uses the equations of |Gimenez| ( |2006[ ). We make this as- 
sumption here because our modeled subsample includes 
only well-detached binary systems. 

The derived parameters for a subset of 18 binary sys- 
tems are presented in Table [l] 3 We have analyzed 
the 18 most promising detached syst ems where a RG 
likely belongs to the EB (see Section 4.1). Thirteen of 
them present RG oscillations, and the other five are the 
longest-period systems where no detectable pulsations 
are observed. This subset actually corresponds to the 
systems with the longest periods and with no ide ntified 
contamination from nearby stars (see Section [3^. The 
light curve data and models are shown in Figures |14| and 

ma 

3.4.2. Search for eclipse timing variations 
As discussed in Section [231 a significant number of the 



candidate RG/EBs are probably not bona-fide RGs in bi- 
nary systems, but are more likely the result of blendings 
from a nearly aligned RG and EB. A further possibil- 
ity is the presence of RGs in triple (multiple) systems in 
which t wo "small" stars mutually eclipse a prim ary com- 
ponent QDerekas et al.||2011| [Carter et al.||2011[ ). Within 
the triple-system hypothesis, it the RGs orbit is elliptical 
and the distance of the RG to the pair of eclipsing stars 

3 An online sortable version of this table is available at 
http:/ /nsol2. nmsu.edu/solarstorm/index.php 
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Figure 14. Modeled light curves for 13 of the 14 longest-period 
RG/EB pulsating systems. In each case, the secondary eclipse 
(defined here as the RG eclipsing the secondary star) has been 
set to an orbital phase of zero. The primary eclipses (secondary 
star eclipsing the RG) have also been aligned, and the given b 
value indicates the phase of the primary eclipse with respect to the 
secondary. (A EB with a circular orbit would have b = 0.5.) The 
y-axis for each panel is the normalized relative flux. The best-fit 
parameters of the models are given in Table ^ 
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Figure 15. Modeled light curves for the five longest-period 
RG/EBs with no RG pulsations. As in Figure [l4| each secondary 
eclipse (defined here as the RG eclipsing the secondary star) has 
been set to an orbital phase of zero. The primary eclipses (sec- 
ondary star eclipsing the RG) have also been aligned, and the given 
b value indicates the phase of the primary eclipse with respect to 
the secondary. The y-axis for each panel is the normalized relative 
flux, and the best-fit parameters of the models are given in Table^ 



is short enough to let tidal forces contribute to the evolu- 
tion of the eclipsing system, we might detect eclipse tim- 
ing variations (ETVs) . This further requires that the RG 
orbital period is, at maximum, the same order of mag- 
nitude as the observation length. In addition, we may 
expect some variations in the eclipse depth and width 
if the precession of the EB system is rapid and strong 
enough to perturb the inclination angle of the EB orbital 
plane. The detection of ETVs is a way to claim a sys- 
tem is at least triple, but the absence of a ETV detection 
cannot be considered proof that a system is not multiple 
since variations may be undetectable or on time scales 
much larger than the observation length. In addition, if 
the RG orbit is almost circular, no ETVs are expected. 
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The presence of ETVs in this particular sample suggests 
that the RG is the third body that perturbs the eclipse 
timing, but this can only be proven with mass determi- 
nations from detailed ETV models and/or the addition 
of radial velocities to the light curve analysis. We note 
that low-mass triple systems are com mon, particularly 



range Kepler's third law to express the stellar mass ratio 
q = M2/M1 as a function of the orbital period P, the 
semi-major axis a, and the mass of the primary star M\\ 



for sy stems with a close-in binary (e.g., |Tokovinin et al. 
2006), so a chance alignment of a triple system with a 
background RG is also possible. 

Modeling ETVs is a common approach to derive 
masses and orbital parameter s for exoplanetary systems 
with at least two planets (e.g.JAgol et al. 1 20051 |Ford et al.| 
2012]|Fabrycky et al.|201^|Stetfen et al.|2013| ) Modeling 
ETVs ol multiple stellar systems is beyond the scope of 
this paper; rather, we first determine whether ETVs are 
detectable and subsequently measure their amplitudes 
and periods. We use one of two methods to measure 
eclipse timings depending on the type of binary system. 
For D and SD sys tems, we fit the p r imary and secondary 
eclipses with the Mandel & Agol (2002) function used 
for fitting exoplanetary transits m the small-planet ap- 
proximation, which acceptably reproduces the behavior 
of stellar eclipses. However, exoplanetary transit func- 
tions are not suitable for OC systems. Therefore, for the 
OC binaries, we measure the minimum intensity time 
by interpolating around the minimum of both eclipses 
with spline interpolation. Such an approach would not 
work for D systems with long orbital periods because 
the eclipse shape is usually strongly modulated by stel- 
lar spots, but it is suitable for SD systems. Once ETVs 
are detected, despite their typically asymmetric profiles, 
we estimate their periods by fitting a simple sine curve 
via the least-square method. 

Eleven systems exhibit ETVs at various signal-to-noise 
ratio levels, as shown in Figure [16] KIC 7990843, 
7031714, 4732015, 10991989, and 7955301 show ETVs 
that are unambiguously coherent between primary and 
secondary eclipses, and for which periods are shorter 
than or near to the total observation length. |Slawson 
et al.| fl201lD note that KIC 7955301 (shown in Figure 
8l presents eclipse depth variations, but based only on 
1-Q2 data they could not determine whether these ef- 
fects were real or due to aperture jitter from quarter to 
quarter. With a longer light curve now available, we see 
this system presents the highest ETV amplitude of about 
four hours. Precession effects are detectable in this sys- 
tem through the shift between primary and secondary 
eclipse timings and the long time-scale evolution of the 
eclipse depths. We also note the unexpected ETVs de- 
tected for contact systems 11135978 and 9181877, which 
show variations on time scales longer than the observa- 
tion length (1400 and 5500 days, respectively) and with 
a mean amplitude of only about 0.83 minutes. This is 
much lower than the 29-minute observing cadence and 
is about half of the ETV standard deviation about the 
fitted sine curve of 1.68 minutes. 

4. DISCUSSION 

4.1. Most RGs do not belong to their associated EB 

4.1.1. Testing binarity with Kepler's law 

We have shown that contamination in the Kepler tar- 
get pixel files confirmed seven systems as false positive 
RG/EBs. To uncover more false positives, we can rear- 
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1. 



(3) 



Note that our convention defines star 1 to be the RG. We 
use RG mas ses a nd radii from the asteroseismic analysis 
(see Section 3.3) with a range of radii for the secondary 
star, 0.1 Rq < R2 < 10 i? , and also consider the rela- 
tive stejj^£^adii_( : R 1+^2) /o> and P from eclipse fitting 
and |Slawson et ah] ( |2011| ). With all of this information 
taken together, we estimate the q required to satisfy Ke- 
pler's third law for the pulsating sample. 

Figure [17| shows the distribution of q in the phase space 
where a 3 /Mi and P are the independent variables. The 
symbols are plotted for the case where R2 = 1R@. The 
horizontal extent of possible values of a 3 /Mi is shown 
over the radius range, where larger radii shift the val- 
ues to the right along the dotted lines. Upon consider- 
ing modest values of i?2, it appears that only 12 of the 
pulsating RGs in candidate RG/EBs may truly belong 
to EBs (not including the ETV systems), as these have 
q < 10. Cases in which the RG would need to be part of 
a binary with a companion star > 10 times more massive 
are considered unrealistic. In this analysis, several sys- 
tems take on a value of q < 0, which is of course unphys- 
ical: this is likely the result of inaccurate estimations of 
the relative stellar radii. The likely 13th candidate (KIC 
5640750) system is not yet conclusive as its orbital pe- 
riod is longer than the observation time. The likely 14th 
candidate (KIC 8095275) does not allow us to estimate 
q since it is a "heartbeat system" for which its orbital 
parameters are unknown and would require specific light 
curve modelin g a nd radial velocities to determine (see 



Section 4.2.2 and lThompson et al. 2012[ ). However, as 



such "heartbeat" light curve features are unique, we are 
confident that it is a bona-fide RG/EB system. 



In summary, the analysis here and in Section 3.1 pro- 
vides strong evidence that only 13 of the 47 pulsating sys- 
tems are true RG/EBs in detached configurations, and 
one is a RG in a non-eclipsing binary system. All of 
these have orbital periods greater than 19 days. All of 
the shorter-period OC and SD systems are either contam- 
inated or the RG is part of a multiple-star system. The 
rest of the D cases must be where the EB and the RG fall 
in the same line-of-sight but may be weakly gravitation- 
ally bound. Finally, we note that testing the likeliness 
of RG/EB candidates is not possible when no RG pulsa- 
tions are detectable. However, true RG/EBs may belong 
to this sample. In particular, the five non-pulsating RGs 
with the longest orbital periods (shown in Figure 15) are 
almost certainly true RG/EBs, as they have light curves 
very similar in appearance to the 13 pulsating bona-fide 
RG/EBs. 

4.1.2. Candidate multiple- star (> 2) systems 

We find eleven stars with ETVs that suggest the RG 
is part of a multiple system composed of a close-in EB 
and a more distant RG in an elliptical orbit. This con- 
figuration, a "hierarchical triple system with t wo low- 
mass stars," was also detected for two c ases in |Carter 



et al. (2011) and |Derekas et aL (2011). In Derekas et al. 
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Figure 16. Eclipse timing variations (ETVs) in hours of the 12 systems that show detectable variations. The y-label "O-C" stands for 
(observed — computed) eclipse epoch. The timings are plotted for the deepest eclipses (black plusses) and for the shallow eclipses when 
detectable (grey dots). Red lines indicate sine-curve fitting on the deepest eclipses. For KICs 3532985, 8255058, and 6762188, the sub-boxes 
zoom in for a more detailed view of the ETVs. 

the disk of the primary component (a subgiant). No 
solar-like oscillations were observed in either case, and 
so here we report the first time global p-mode oscilla- 
tions are detected for an RG in such a system. In the 
eleven ETV detections, only the close-in eclipses are ob- 
served and the RG's presence is indirectly deduced as the 
perturber. We find this in seven D, two SD, and two OC 
cases. The system KIC 4758368 is the only out of the 
eleven where we do not detect any RG oscillations (see 
Table 0. 

4.2. EB and multiple- star candidates 

We classify the RGs displaying oscillations and pre- 
sumably belonging to a double or triple system into one 
of three categories. First, fundamental cases are those 
that deserve to be studied in more detail in the future, 
to be used as cornerstones for testing stellar-evolution 
models, and for which current Kepler data are sufficient 
for accurate eclipse and asteroseismic modelings. Sec- 
ond, the promising cases are those for which eclipse or 
asteroseismic modeling require additional Kepler obser- 
vations to be considered fundamental. Third, intriguing 
cases are those in which we cannot rule out the detected 
signal as a false positive, and where more data could lead 
to unexpected discoveries. 

4.2.1. Fundamental cases 

The fundamental cases are EBs and one hierarchical 
triple system, with orbital periods over 19 days. Due to 
the high signal-to- noise ratio of their oscillation patterns, 
these systems are suitable for precise modeling of their 




Figure 17. Mass ratio q as a function of orbital period and a 3 /Mi 
for the pulsatin g RG/EB candida te sample based on the orbital pa- 
rameters in |Slawson et aL] ( |2011| > and those obtained in this work 
from asteroseismology. tor all systems, the symbols are plotted as- 
suming R2 = Rq, and the dashed lines show the extent of a 3 /Mi if 
the secondary star's radius spans 0.1 Rq < R2 < 10 Rq. The color 
at each system's position represents the mass ratio q. Contours of 
constant q = 10 and 100 are shown by the diagonal lines. The size 
of each symbol is proportional to the pulsating RG stellar radius, 
Ri , measured asteroseismically. The pulsating OC systems are not 
included here. The colors reflect the huge variation of mass ratio 
across the sample, and clearly distinguish which systems are likely 
to be a bona-fide RG in an EB (yellow, i.e., q ~ 1) from those that 
are unlikely (blue, i.e., q > 10). 



(2011) the low-mass stars were co-planar with the pri- 
mary c omponent (a RG) and all eclipses were visible, 
while in |Carter et al.| ( |2Qll| ), each star separately eclipses 
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interior properties. The candidates with highly eccentric 
orbits are interesting for studying tidal interactions in 
multiple-star systems. Their modeled light curves are 
shown in Figure [14] We list them here, sorted by the 
spectral type of the companion star: 

M or K dwarf companion. We identify two cases 
where the EB is composed of a RG and a smaller, 
cooler companion, rendering the primary eclipse 
deeper than the secondary eclipse. 

- KIC 8702921 is a 19-day EB, which shows 
clear RG oscillations typical of the RGB. Its as- 
teroseismic radius and mass are 5.1 ± 0.1 Rq and 
1.4 ± 0.1 M©. We find that the secondary star is 
compatible with an M dwarf, as its mass is esti- 
mated as 0.6 ± 0.3 M© with radius 0.44 ± 0.01 i? , 
and an effective temperature of 2600 K. We see 
variations in the light curve that could be due to M- 
star activity (i.e., spots), tidal distortions, and el- 
lipsoidal and Doppler beaming. We observe a 97.8 
day (almost precisely five times the orbital period) 
modulation in the light curve of rather large am- 
plitude, which is comparable to the eclipse depth. 

- KIC 5308778 has a 41-day orbital period and 
low signal-to- noise ratio RG oscillations, corre- 
sponding to an asteroseismic radius and mass of 
9.4 ± 0.8 R Q and 1.3 ± 0.3 M©. From eclipse 
modeling, the companion's radius and mass R^ — 
0.6 ± 0.1 R Q and M 2 = 0.3±oI M Q ^ with an M 
dwarf, while its effective temperature T 2 = 4140 K 
fits with a K dwarf. Its light curve presents photo- 
metric variations of amplitude up to 6%, which is 
17 times larger than the deepest eclipse. These 
variations are quasi-periodic with periodicity al- 
most equal to the orbital period. The absence of 
any orbital eccentricity suggests that the photo- 
metric variations are due to features on the RG, 
in a system that is tidally locked on a circularized 
orbit. The periodic fluctuations appears incompat- 
ible with spots on the companion star since its rel- 
ative brightness is 0.2 % of the RG. 

G or late F dwarf companion. Three systems are 
characterized by a RG with a G or late type F 
star on an eccentric orbit with typical e ~ 0.2, and 
light curves that present large variations. 

- KIC 8430105 is a 63-day period EB with clear 
RG oscillations that correspond to an asteroseismic 
radius and mass of 7.4 ± 0.3 R@ and 1.2 ± 0.1 M©. 
No mixed modes are identified, but the large sep- 
aration value and the seismic mass estimate sug- 
gest this RG cannot belong to the red clump and 
is more likely an RGB star. From eclipse model- 
ing, the orbit appears to be eccentric as e = 0.26, 
and we estimate the companion's radius, mass, and 
effective temperature to be R2 = 0.83 ± 0.04 

M 2 = 0.9 ± 0.4 M©, and T 2 = 5960 K. Thus, we 
likely have an RG in orbit with a solar analog. A 
strong variation in the light curve of (almost pre- 
cisely) two times the orbital period and five times 
the eclipse depth is found. At the bottom of the 
primary eclipse (the G star passing in front of the 
RG), we observe sharp peaks that can reach up to 
half of the eclipse depth and last about of third of 
the eclipse time. These could be strong flares on 



the solar analog, which seem unlikely, or hot spots. 

- KIC 10001167 is a 120-day EB with clear RG 
oscillations that are characteristic of an RGB star 
of a 14.0 ± 0.8 R Q and 1.1 ± 0.2 M©. Modeling 
the eclipses indicates the orbit is eccentric with 
e = 0.16, and the companion is likely a GO or 
F9 star with radius R2 = 1.1 ± 0.1 i?©, mass 
M 2 = 1.0 ± 0.6 M , and effective temperature 
T 2 = 6090 K. 

- KIC 9970396 presents a 235-day period and 
clear modes of a RG with 8.3 ± 0.3 Rq and 1.3 ± 
0.1 M . The orbit appears to be eccentric with 
e = 0.20, and its companion is likely a late F star of 
radius R 2 = l.liO.l i? , mass M 2 = 1.2±0.2 M©, 
and effective temperature T 2 = 6060 K. In contrast 
to the two previous cases, its light curve does not 
exhibit high photometric variations. 

F type companion. These four bona-fide RG/EBs are 
likely composed of an RGB star with an F type 
companion on an eccentric orbit (e = 0.23 — 0.67), 
with orbital periods ranging from 175 to 408 days. 
In these cases, the secondary F star's mass is es- 
timated to be nearly these same or even slightly 
larger than the RG star mass. However, the er- 
ror bars on the masses are large so that the pa- 
rameters are consistent with the RG mass being 
larger than the secondary star mass, as would be 
expected if the stars formed together and the more 
massive RG finished its main-sequence phase first. 
However, it is also possible that the RG could lose 
several tenths of a solar mass in a wind during the 
RG phase, or could have transferred mass to the 
companion since the orbits are highly eccentric, so 
that at the present time the RG has a slightly lower 
mass than its companion. 

- KIC 9540226 presents a 175.5-day orbital pe- 
riod and clear RG oscillations that are character- 
istic of an AGB or RGB star with 14.0 ± 0.7 R Q 
and 1.6 ± 0.2 M©. Eclipse modeling indicates the 
orbit is eccentric with e = 0.39, and the compan- 
ion is likely an F star of radius R2 = 1.4 ± 0.1 i?©, 
mass M 2 = 1.5 ± 0.7 M©, and surface temperature 
T 2 = 6920 K. 

- KIC 5786154 presents a 197.9-day period and 
shows clear oscillations of a a RG of 12.7 ± 0.6 Rq 
and 1.4 ± 0.2 M©. The orbit appears to be eccen- 
tric with e = 0.38, and the companion is likely 
an F star of radius R2 = 1.9 ± 0.1 i?©, mass 
M 2 = 1.8 ± 0.7 M©, and surface temperature 
T 2 = 6545 K. 

- KIC 7037405 presents a 207-day period and 
clear modes of an RGB star with 15.0 ± 0.9 Rq 
and 1.4 ± 0.2 M©. The orbit appears to be eccen- 
tric with e = 0.23, and the companion is likely 
an F star of radius R2 = 1.9 ± 0.1 i?©, mass 
M 2 = 1.5 ± 0.7 M©, and surface temperature 
T 2 = 6400 K. 

- KIC 8410637: This system was the only EB 
previously k nown to host a RG displaying global 
oscillations (Hekker et al. 2010). Its light curve 
is one of the most challenging among the set of 
EBs to model because of a high eccentricity (e = 
0.67), coupled with relatively small stellar radii 
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compared to the semi-major axis, (Ri + R 2 )/a = 
3.67 %. We confirm the system is composed of a RG 
and an F star, whose respective sizes and masses 
are (R u Mi) = (11.0 ± 0.5#©, 1.6 ± O.2M ) and 
(R 2 , M 2 ) = (1.6 ± O.LR©, 1.8 ± O.7M ). 

Probable £-Scuti companion. The system KIC 
4663185 presents a 57-day orbital period and 
is the only from our sample that has a power 
spectrum with two clear sets of oscillations. On 
the one hand, RG oscillations are observed at 
^max = 23 /iHz: this corresponds to a rather 
massive RGB star with radius 18.2 ± 1.2 R@ and 
mass 2.2 ± 0.4 M . On the other hand, a second 
oscillation spectrum is observed at ^ max ~ 132 /iHz 
(Fig. 18). These mode amplitudes and widths 
are too nigh and narrow, i.e., the lifetimes are too 
long to match with a RG oscillation spectrum. 
The oscillation spectrum with ^ max — 132 /iHz is 
consistent with a J-Scuti or /3-Cep star p-mode 
spectrum. However, a /3-Cep star is probably too 
massive (~ 10 M ) to have evolved together with 
an RGB star of current mass 2.2 M . Regarding 
the light curve, only one type of eclipse is detected 
and they appear to be shallow and grazing. We 
associate the detected eclipses with the secondary 
transit, i.e., the companion star is eclipsed, since 
the 5-Scuti is likely to present a higher surface 
brightness than the RG. With this hypothesis, 
we put an upper on the companion's size of 
R 2 < 1.5 ±0.1 R®. The companion's size estimate 
is consistent with a main-sequence star of mass 
~ 1.6M , which fits with the interpretation that 
the secondary could be a £-Scuti star. In the 
future, we will obtain ground-based spectroscopic 
observations, and attempt to match the oscillation 
frequency spectrum with stellar pulsation models 
to confirm this hypothesis and better constrain 
the parameters for the secondary. 

RG companion. The system KIC 9246715 is a 171.3- 
day EB showing RG oscillations with a low signal- 
to-noise ratio compared to the other oscillating 
RGs, despite it being the brightest system (K p — 
9.266). From eclipse modeling it appears that this 
system has an eccentric orbit (e = 0.35) and is 
composed of a pair of RGs of similar radii and 
temperature (R2/R1 = 0.77, T 2 /T 1 = 1.03). How- 
ever, we detect only one oscillation pattern (^ max = 
102 /iHz), whose asteroseismic parameters corre- 
spond to a RG of 7.7 ± 0.4 R Q and 1.7 ± 0.3 M . 
The large separation indicates it is a RGB star. 
From this, we deduce the companion's radius and 
mass to be R 2 = 5.3±0.3 R Q ad M 2 = 0.8±0.7 R Q . 
No pulsations from the companion star are de- 
tected: it likely pulsates with a lower signal-to- 
noise ratio in the same frequency range, since the 
companion's frequency at maximum amplitude is 
expected to be ^ max = 96^g4 3 /iHz. In addition, the 
autocorrelation of the light curve reveals the pres- 
ence of a sine modulation with about a 3-day period 
with a mean peak-to-peak amplitude of about 100 
ppm. 

Hierachical triple system. The system KIC 7955301 
displays high signal-to-noise ratio ETVs (4-h am- 



plitude, Figure 16) and clear RG pulsations, typical 



of an RGB star, with 1.2±0.1 M , 5.9±0.2 R Q , and 
a core rotation period of 30 days. From the eclipse 
timings, we are able to infer that the RG takes 210 
days to orbit the more compact 15-day period EB. 
The components of this system are likely to inter- 
act strongly as indicated by the high and complex 
amplitude of their ETVs. The asymmetrical shape 
of the ETV curve with respect to a sine curve sug- 
gests the RG orbit is highly eccentric. In addition, 
the variable eclipse depth as function of time, on 
a period of about 1800 days, suggests that the or- 
bital plane of the pair of EB is precessing due to 
tidal interactions with the RG. 

4.2.2. Promising cases 

We sort these potentially fundamental cases by increas- 
ing orbital period. 

KIC 5179609: This system has a value of q in a plau- 
sible range (0.6) with a 43.9-day period. Unfor- 
tunately, the global modes of the RG are largely 
at frequencies higher than the Nyquist frequency, 
so we are only able to properly measure the mean 
large separation. Asteroseismic scaling laws infer 
the RG's mass and radius to be 1.3 ± 0.1 M and 
3.7 ±0.1 Rq. In addition, we do not find a de- 
tectable signature of the secondary eclipse in this 
EB. Given that the standard deviation of the light 
curve is 164 ppm and the primary eclipse depth is 
rather significant (0.92 %), the ratio R 2 T% / (R\Tf-\- 
R 2 T$) is lower than 1.78%. This is a result of a 
small or cool companion star. For simplicity, we 
can model the system as an exoplanetary system, 
i.e., as a dark planet eclipsing a star. With such 
an assumption, the radius of the companion would 
be O.4i? , which could then be an M, L, or brown 
dwarf. One quarter or Kepler data at short cadence 
would allow us to better constrain this system and 
make it a cornerstone of for testing asteroseismol- 
ogy. It deserves further study as it hosts the small- 
est companion star of our sample. 

KIC 8095275: Highly eccentric but non-eclipsing bi- 
nary systems have been observed by Kepler, and 
are commonly referred to as "heartbeat stars" due 
to the resemblance of th eir light curves to an elec- 
trocardioeram (Fig. [18) I Welsh et al.|2011| [Thomp- 
son et al. 2012. In these cases, the photometric 
variability is due to tidally induced distortions gen- 
erated by the companion star as it passes close to 
the pericenter. The system KIC 8095275 presents 
a 23-day orbital period, a clear heartbeat signal 
of about 1 % amplitude, and RG oscillations that 
correspond to a radius of 7.6 ± 0.2 R@ and mass of 
1.2 ± 0.1 Mq. A more detailed model of the light 
curve coupled with radial velocity measurements 
should aid us in characterizing this system. We 
note that it was also detected in parallel to this 
work and publicly reported on the Internet, even 
though it is not yet mentioned in any peer-reviewed 
paper. 4 

4 http://keplerlightcurves.blogspot.com/2012/09/three-giant- 
heartbeats.html 
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KIC 4732015 and 10991989: These D systems both 
present clear ETVs with periods on the order of 
the observation duration (1050 and 544 days, re- 
spectively) with a rather small amplitude (6.9 and 
3.9 min). Both systems are probably hierarchical 
triple systems that have a more distant RG orbiting 
a pair of close main-sequence stars (0.93 and 0.97 
day periods, respectively). However, these systems 
are significantly different physically. KIC 4732015 
is an RGB star and is the biggest RG in our sample 
(#i,Mi) = (33.0 ±5.1 i? ,1.5 ±0.5 M©), while 
KIC 10991989 is a quiet, massive RC2 star with 
(R u Mi) = (9.6 ± 0.5 #©,2.5 ± 0.4 M©). 

KIC 5640750: This system is the only one whose or- 
bital period is longer than the total observation 
time (Q0-Q13; P > 1149.5 days). It appears that 
the eclipse shape is similar to those of RGs orbit- 
ing F type companions. Modeling the light curve 
yields a degeneracy between the orbital period and 
the eccentricity because we see only one primary 
and one secondary eclipse so far. For our pur- 
pose, however, the temperature and radius ratios 
may still be estimated even if the mass and semi- 
major axis are not yet. To model the light curve, 
we assumed the orbital period to be twice the time 
between the primary and secondary eclipses and 
set the argument of the periastron to uj = 7r/2, 
because such an assumption leads to the lower 
boundary of the eccentricity. From asteroseismic 
measurements, the RG has radius and mass of 
CRi,Mi) = (14.3 ±0.8 i? ,1.4 ±0.2 M©). Cou- 
pling this result with the preliminary eclipse mod- 
eling indicates that the companion seems to be an 
F star with radius R2 = 1.8 ± 0.1 Rq and tem- 
perature T 2 = 6497 K, on an elliptical orbit with 
e > 0.18. The high value of the RGs asteroseismic 
mass (3.2 ±1.0 M©) suggests its period and semi- 
major axis are largely underestimated, which also 
suggests a higher eccentricity. 

4.2.3. Intriguing cases 

KIC 3532985 and 6762188: These D systems both 
present ETVs with short periods (29 and 35 days, 
respectively) of small amplitude (3.2 and 6.5 min, 
respectively), with asymmetrical shapes. At first 
glance, the ETVs look like artifacts, but several 
methods to determine the eclipse timing were used 
and all yielded ETVs in these systems but not in 
the other targets. The RGs associated with these 
systems are red clump stars: KIC 3532985 is a RC2 
with (R u Mi) = (7.8±0.2 #©, 1.7±0.1 M©), while 
KIC 6762188 is a solar-mass RC1 with (R u Mi) = 
(10.3 ±0.5 # ,1.O ±0.1 M ). 

KIC 8255058: This D system presents the most puz- 
zling ETV features from our sample. It appears 
similar to KICs 3532985 and 6762188, but its ETVs 
actually split into two components. We observe 
a rapid oscillation of period equal to the EB or- 
bital period, with 5-minute amplitude, as well as a 
long period (868 days) of similar amplitude (6.7 
min). The RG is an RGB star and is one of 
the smallest RGs in our sample with (i?i,Mi) = 
(3.9 ±0.1 i? ,1.4 ±0.1 M ). 



.x 10 



o 1 — 

100 



,,11. 




120 130 

v (/xHz) 



140 



150 



160 




0.992 



10 15 
Time (days) 



Figure 18. Top: oscillation spectrum of the RG's companion of 
KIC 4663185, which is likely a (5-Scuti star. Bottom: folded light 
curve of the "heartbeat" star KIC 8095275, which is the only from 
our sample to be a RG in a non-eclipsing binary. Black crosses 
indicate the folded signal, while the red line indicates the same 
rebinned in 30-min bins. 

KIC 7690843 and 7031714: These SD systems both 
present clear ETVs with 75 and > 1000 day peri- 
ods, respectively. The amplitude of the ETVs are 
weak and on the order of a few minutes, i.e., shorter 
than Kepler's long-cadence sampling. However, 
KIC 7031714 presents evidence for stellar crowd- 
ing and blending. 

KIC 11135978 and 9181877: These OC systems dis- 
play very slowly- varying ETV trends (> 1400 day) 
with low signal-to- noise ratios in amplitude. One 
would not expect that a third body could signif- 
icantly perturb the orbital equilibrium of the two 
contact stars. 

KIC 7377422, 4569590, 3955867, 9291629, & 7943602: 

These detached systems present orbital periods, 
eclipse depths, and eclipse shapes similar to the 
eight strongest RG/EB candidates where we detect 
RG pulsations (see Figures 14 and 15). Oddly, 
the light curves of these four systems ao not show 
pulsations. Their peculiarity, with respect to the 
pulsating RG/EBs candidates, consists of high 
quasi-periodic photometric variability (between 10 
and 30%). This suggests that either their modes 
are buried in the noise of photometric variations, 
or that their mode frequencies are above the 
Nyquist frequency. Alternatively, the fact that no 
solar-like oscillations are detected in these objects 
that display signific a nt stellar variability agrees 
with |Chaplin et al.| ( |2011[ ), which suggests that 
stellar activity inhibits the amplitude of solar-like 
oscillations. This scenario does assume that the 
stellar variability originates from the RG and not 
the companion, which is unknown. 
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5. CONCLUSION AND PROSPECTS 

We have identified a set of potentially useful targets 
to test theories of stellar evolution. We are confident 
that 13 (of 70) systems cross-listed from the RG and 
EB Kepler catalogs are bona-fide pulsating red giants in 
eclipsing binaries. One is a red giant in a non-eclipsing 
binary system, and an additional five are likely red giants 
in eclipsing binaries even though we do not detect their 
pulsations. It is likely that at least 11 other systems 
are candidates for belonging to three-body configurations 
composed of a pair of eclipsing main-sequence stars and 
a red giant. 

Oddly, we detect no oscillations in 23 of the red gi- 
ant stars in candidate RG/EB systems across all classes 
of binaries. We identify several reasons that may explain 
this. First, the star could be misidentified and is not a gi- 
ant. Second, the treatment used for removing the eclipse 
signature may suppress any oscillations at very low fre- 
quency. Also, the giants could actually be subgiants or 
low-luminosity red giants pulsating at frequencies be- 
yond the Nyquist frequency. Finally, some phenomena 
in the interactions of a multiple system may yield a large 
damping e ffect of the global oscillations in the RG ( | Fuller | 
et al.||2Q13 ). Therefore, in total, we believe many oi the 
RG/EB candidates deserve to be followed-up with short- 
cadence Kepler observations. This would allow global 
modes to be detected in a red giant for which ^ max is 
larger than the long-cadence Nyquist frequency. Ideally, 
this would allow for pulsations of the main-sequence com- 
panions to be measured as well; however, the contribu- 
tion to the total light from the companion star is only 
several percent in the best cases, making such detections 
very challenging in practice. 

Spectroscopic measurements from the ground will cer- 
tainly help in understanding these systems too. First, 
the identification of overlapping spectra could indicate 
whether the system is indeed composed of a RG and a 
main-sequence star. Second, the measurement of radial 
velocities is a way of extracting accurate masses for the 
system that can be cross-checked with the asteroseismic 
inferences. Third, it is a way to refine the estimate of 
the stellar parameters which are currently based on color 
photometry. We have already started observations of a 
subset of these 70 stars with the ARCES echelle spec- 
trometer (resolution R = 30 000) at the Apache Point 
Observatory (APO), New Mexico. Other targets from 
our sample are currently being observed by the APOGEE 
spectrometer on the SLOAN Digital Sky Survey (SDSS) 
telescope at APO, in the context of the APOKASC pro- 
gram to support Kepler observations in asteroseismology. 
This effort is coordinated by teams from the APOGEE 
project and the Kepler Asteroseismic Science Consor- 
tium (KASC). In addition, the KASC Working Group 8 
is studying in detail the systems KIC 8410637, 5640750, 
and 9540226. If more precise estimates of the stellar 
parameters can be obtained and coupled with eclipse in- 
formation and asteroseismology, these RG/EB systems 
have the potential to become some of the most accurately 
studied stars. 
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Table 2 

Derived properties of the 47 pulsating systems from combined asteroseismic 
and light-curve analysis. 



Result | KIC parameters | Asteroseismic analysis and combined output 
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Note. - Only the first 13 EB lig ht curves were modeled, as these are the most likely true RG-EB systems. See also 

http: / /nsol2. nmsu.edu/solarstorm/index.php 

a Core rotation period determined from mixed-mode analysis 

k Asteroseismic RG classification: Red-giant branch (RGB); Asymptotic red-giant branch (AGB); Red clump 1 (RC1); Red clump 2 (RC2) 



